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The numerical technique of discrete light-cone quantization (DLCQ) is applied to a 
single-fcrmion truncation of Yukawa theory in four dimensions. The truncated theory 
is regulated by three Pauli-Villars bosons, which are introduced directly in the DLCQ 
Fock-state basis. A special form of the Lanczos diagonalization algorithm is used to 
handle the indefinite metric. Renormalization is done nonperturbatively. 

1. Introduction 

Methods for the nonperturbative numerical solution of light-cone-quantized quan- 
tum field theories have progressed to the point where they are applicable in four 
dimensions. The use of light-cone coordinates™ makes possible a meaningful Fock- 
state expansion, in which no disconnected vacuum contributions appear. The Fock- 
state wave functions in the expansions are obtained by solving a Hamiltonian eigen- 
value problem.0 A standard technique for solving such a problem is discrete light- 
cone quantization (DLCQ)BB The wave functions are evaluated at discrete momen- 
tum values p + = E + p z = nit/L, p± = (p x ,p y ) = ri±ir/Lj_, with n, n x , and n y 
integers and L and L± length scales. The coupled integral equations that comprise 
the eigenvalue problem become a matrix diagonalization problem where trapezoidal 
sums approximate integrals. 

To properly formulate such a matrix problem as an approximation to a four- 
dimensional field theory, one must include a regularization scheme and perhaps 
additional cutoffs that yield a finite matrix problem. One must also include a 
renormalization scheme to determine bare parameters. These steps have been car- 
ried out for simple modelslffl and are now being applied to Yukawa theory!! The 
regularization scheme is based on the introduction of Pauli-Villars particles ,0 in- 
cluding some with negative norm, and of a simple mass counterterm. The theory 
is then finite before discretization. A cutoff on the light-cone energy (to 2 J rp\)/p + 
is used to limit the transverse momentum range and produce a matrix representa- 
tion of finite size. Renormalization is done nonperturbatively, by fixing computable 
quantities to "data." 

The original theory is recovered in the following sequence of limits. First, the 
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numerical limit of infinite longitudinal and transverse resolutions is taken at fixed 
cutoff and fixed Pauli-Villars masses. Next, the cutoff is removed, and, finally, the 
Pauli-Villars masses are taken to infinity. Calculations done to datetffl do not show 
a great sensitivity to numerical resolution, above a modest threshold, so that the 
basis sizes used in the matrix problem have been manageable. The largest basis 
used is approximately f0.5 million states. 

2. Yukawa Theory 

The light-cone Hamiltonian for Yukawa theory is given by McCartor and Robertson.i 
We work with a single-fermion truncations in which no pair creation or annihila- 
tion terms appear and for which an eigensolution is sought only in thcxine-fermion 
sector. An analysis of the one-loop fermion self-energy then showsBcl that three 
Pauli-Villars bosons are necessary and that their couplings to the fermion are fixed 
as functions of their masses by three algebraic conditions. Two of these bosons 
must be negatively normed. A mass counterterm is also included. The singularity 
in the instantaneous fermion contribution is canceled by_ addition of an effective 
interaction patterned after the contribution of a Z graph.Q 

Renormalization is done by holding fixed the mass of the dressed fermion state 
and by fixing the value of the expectation value for :0 2 (O):, where <fi is the boson 
field operator. The eigenvalue problem is rearranged so that the coefficient of the 
mass counterterm becomes the eigenvalue. In this form the eigenvalue problem is 
solved simultaneously with the condition on (:</> 2 (0):) by iterating in the value of 
the bare coupling. This determines the bare mass and bare coupling as functions of 
the numerical parameters and the regularization parameters. The Fock-state wave 
functions are also obtained. 

3. Numerical Methods and Results 

The matrix eigenvalue problem is solved with a variant of the biorthogonal Lanc- 
zos methodS designed specifically for an indefinite metric! This iterative method 
requires stopping criteria, which we take to be convergence of the eigenvalue and of 
parts of the boson-fermion wave function, as well. Because the method generates 
several eigenvalues, we also need criteria for selecting the state of interest. This 
state is the ground state, but due to the indefinite metric, it is not necessarily the 
state of lowest mass. The criteria used for selection include the following: a positive 
norm, a real eigenvalue, absence of nodes in the parallcl-hclicity boson-fermion wave 
function, and a relatively large bare-fermion probability. A starting point for the 
iterations is generated with use of high-order Brillouin-Wigner perturbation theory. 

An important check on the calculation is found in the antiparallel boson-fermion 
wave function, which due to J z conservation must be in an L z — 1 state. The 
calculation does not assume this symmetry but instead computes the wave function 
at all n x and n y values. This wave function is found to have the correct symmetry. 

Given a method for the computation of Fock-state wave functions, any number of 
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interesting quantities can be subsequently computed. For example, matrix elements 
of the fermion current operator yield form factors for the dressed state.Ell 

4. Future Work 

The techniques described here are applicable to quantum electrodynamics and pos- 
sibly quantum chromodynamics. JLn the latter case one would need to use a formula- 
tion such as that of Paston et alS-n where an appropriate number of ghost particles 
and counterterms are inserted, or perhaps a theory with (broken) supersymmetry. 

In the short term, a more complete investigation of Yukawa theory can be made. 
The two-fermion sector of the no-pair version is of interest because one can consider 
true bound states and scattering states. The full theory is also interesting because 
pair terms bring additional divergences and renormalization of the boson mass. 
These cases can be pursued with direct extensions of present methods. 

Acknowledgements 

The work reported here was done in collaboration with S.J. Brodsky and G. McCar- 
tor and was supported in part by the Department of Energy, contract DE-FG02- 
98ER41087, and by grants of computing time from the Minnesota Supercomputing 
Institute. 



References 



1. P.A.M. Dirac, Rev. Mod. Phys. 21 (1949) 392. 

2. For reviews, see S.J. Brodsky, H.-C. Pauli, and S.S. Pinsky, Phys. Rep. 301 (1998) 299; 
J.R. Hiller, to appear in the proceedings of the Tenth International Light-Cone Meeting 
on Nonperturba tive QCD and Hadron Phenomenology, Heidelberg, June 12-17, 2000, 
hep-ph/0007309[ 

3. H.-C. Pauli and S.J. Brodsky, Phys. Rev. D 32 (1985) 1993; 32 (1985) 2001. 

4. S.J. Brodsky, J.R. Hiller, and G. McCartor, Phys. Rev. D 58 (1998) 025005. 

5. S.J. Brodsky, J.R. Hiller, and G. McCartor, Phys. Rev. D 60 (1999) 054506. 

6. J.R. Hiller, to appear in the proceedings of the CSSM Workshop on Light-Cone QCD 



and Nonper turbative Hadron Physics, Adelaide, Australia, December 13-21, 1999, hep- 
ph/0002222} 



7. W. Pauli and F. Villars, Rev. Mod. Phys. 21 (1949) 4334. 

8. G. McCartor and D.G. Robertson, Z. Phys. C 53 (1992) 679. 

9. C. Bouchiat, P. Fayet, and N. Sourlas, Lett. Nuovo Cim. 4 (1972) 9; S.-J. Chang and 
T.-M. Yan, Phys. Rev. D 7 (1973) 1147; M. Burkardt and A. Langnau, Phys. Rev. D 
44 (1991) 1187. 

10. C. Lanczos, J. Res. Nat. Bur. Stand. 45 (1950) 255; J. Cullum and R.A. Willoughby, 
in Large-Scale Eigenvalue Problems, eds. J. Cullum and R.A. Willoughby, Math. 
Stud. 127 (Elsevier, Amsterdam, 1986), p. 193. 

11. S.J. Brodsky and SD. Drell, Phys. Rev. D 22 (1980) 2236; J.R. Hiller, in the pro- 
ceedings of the TJNAF workshop on the Transition from L ow to High Q Form Factors, 
Athens, GA, September 17, 1999, p. 193, |hep-ph/9909471 . 

12. S.A. Paston, V.A. Franke, and E.V. Prokhvatilov, Theor. Math. Phys. 120 (1999) 
1164. 



3 



